Analog of Magnetoelectric Effect in High-T c Granular Superconductors 
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We propose the existence of an electric-field induced nonlinear magnetization in a weakly coupled granular su- 
perconductor due to time-parity violation. As the field increases the induced magnetization passes from para- to 
dia-magnetic behavior. We discuss conditions under which this effect could be experimentally measured in high 
temperature superconductors. 



By applying a constant external electric field through 
an insulating layer, recent experimental studies 
have found an enhancement in the critical current of 
ceramic high-T c superconductors (HTCS). It has been 
shown that for any field-induced effects to be experi- 
mentally measurable, the electrostatic screening length 
should be larger than the superconducting coherence 
length. This happens to be the case in HTCS, which 
have a reduced carrier number density To account 

for the unusual behavior of the field dependent critical 
current observed in HTCS ceramics however, it was 

argued that some other mechanisms (in addition to the 
conventional carrier number density modification) have 
to be taken into consideration. In particular, it was sug- 
gested that an applied electric field may actually induce a 
substantial change into the original weak-links structure 
of granular samples || . 

In this letter we predict yet another interesting effect 
which can be produced by an external electric field ap- 
plied to a granular material. This effect is similar to 
(but physically different from) the magnetoelectric effect 
(MEE) seen in magnetic ferroelectrics (like antiferromag- 
netic BiFeO-i) 0. It has also been discussed in the con- 
text of some exotic nonmagnetic normal metal conduc- 
tors, with the symmetry of a mirror isomer ||, and pyro- 
electric superconductors (where the supercurrent passing 
through a metal of a polar symmetry is assumed to be ac- 
companied by spin polarization of the carriers |^]). The 
effect discussed here entails the electric field generation 
of a magnetic moment (paramagnetic in small to dia- 
magnetic for large electric fields) due to superconducting 
currents that circulate between grains. The basic phys- 
ical reason for the appearance of this effect is that the 
applied electric field induces a magnetization M, that 
changes its sign (M — > — M) under a time-parity trans- 
formation. There have been other microscopic mecha- 
nisms proposed to have time-parity violation effects in 
HTCS. One that attracted significant experimental at- 
tention was a theoretical consequence of the anyon pic- 
ture of HTCS 1 10 1 . Another proposition was based on a 
spin-orbit coupling 0. Up to now, however, the experi- 
mental search for their signatures have not bared fruit. 

In the case of weakly-coupled superconducting grains, 
the phenomenological reason for the time-parity viola- 



tion comes from the fact that the total free energy F for 
the material exhibiting both magnetic and electric prop- 
erties contains the term a mn f mn (E, H), with the coeffi- 
cients a mn 7^ (here {m,n} = x,y,z). In the standard 
MEE, the function f mn = E m H n , which leads to the cor- 
responding linear effect, which is either an electric-field 
induced magnetization M m (E) = dF/dH rn — a mn E n 
(in zero magnetic field) or a magnetic-field induced po- 
larization P m (H) = dF/dE m — a mn H n (in zero elec- 
tric field). Since H (and M) changes its sign under 
time-parity transformation while E (and P) remains 
unchanged, an electric field induced magnetization will 
break time-parity symmetry even in zero magnetic field 
applied. As we show below, in our case the symmetry 
breaking term in F, represented by a nonzero coefficient 
a mn , has a more general nonlinear form for f mn . We 
recall that the standard linear MEE can appear when an 
external electric field E interacts with an inner magnetic 
field Hdm of the Dzyaloshinskii-Moriya (DM) type | flT| . 
The DM interaction leads to a term, apart from the stan- 
dard isotropic term in the Heisenberg Hamiltonian, of the 
form Hdm = J2i j A,j 1 (Si A Sj), where Si is a Heisen- 
berg spin and the constant vector Di j arises from the 
spin-orbit coupling. An analogous situation occurs in 
our case, as we describe below. To see how we can get 
a nonzero a, or equivalently a DM type interaction in 
a granular superconductor, we model a HTCS ceramic 
sample by a random three-dimensional (3D) overdamped 
Josephson junction array. This model has proven to be 
useful in describing the metastable magnetic properties 
of HTCS |f^,|l3|. In thermodynamic equilibrium, this 
model has a Boltzmann factor with a random 3D-XY 
model Hamiltonian. Specifically, the general form of the 
Hamiltonian (describing both DC and AC effects) reads 



(1) 



Ej( r ij) is the 
?i — fj the sep- 
aration between the grains; the gauge invariant phase 
difference is defined as 



Here {i} = ri is a 3D lattice vector; 
Josephson coupling energy, with fi 



Mt) = MO) - Aij(t), 
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where <pij (0) = cf>i — 4>j with <f>i being the phase of 
the superconducting order parameter; Ay(t) is (time- 
dependent, in general) frustration parameter, defined as 



= IT I Ar,t)-dl, 
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with A(f, t) the (space-time dependent) electromagnetic 
vector potential which involves both external fields and 
the electric and magnetic possible self-field effects (see 
below); $o = h/2e is the quantum of flux, with h 
Planck's constant, and e the electronic charge. Ex- 
panding the cosine term, and using trigonometric iden- 
tities we can explicitly rewrite the above Hamiltonian as 
Hxy = - Eij Ej[cos(A itj )Si ■ Sj - sin(A iJ )fc • S, A Sj}, 
where the two-component XY spin vector is defined as 
Si = (cos 4>i, sin^i), and A: is a unit vector along the 



z-axis 1 14 . We see that the second term in this Hamilto- 
nian has the same form as in the DM contribution, and 
thus we can surmise that the time parity will be broken 
by applying an external field (electric or magnetic) to the 
granular system. To bring to the fore this possibility, we 
show below in a simple but yet nontrivial model that this 
is indeed the case. 

One important property of the Hxy Hamiltonian is 
that it is random because Aij itself is a random variable. 
There are different types of Ai.j randomness that can be 
considered Jl3| . For simplicity, in the present paper, we 
consider a long-range (Sherrington-Kirkpatrick-like) in- 
teraction between grains (assuming Ej(jij) = Ej) and 
model the true short-range behavior of a ceramics sample 
through the randomness in the position of the supercon- 
ducting grains in the array (using the exponential dis- 
tribution law P r (rij), see below). Here we restrict our 
consideration to the case of an external electric field only 
but it can be shown that the scenario suggested in this 
paper will also carry through when applying an exter- 
nal magnetic field (which will induce another time-parity 
breaking phenomenon in the granular material; namely, 
magnetic field induced electric polarizability [^5|). Be- 
sides, in what follows we also ignore the role of Coulomb 
interaction effects assuming that the grain's charging en- 
ergy E c Ej (where E c — e 2 /2C, with C the capaci- 
tance of the junction). 

In the case of a granular material, we show here that 
the corresponding time-parity breaking DM internal field 
can be related to the electric field induced magnetic mo- 
ment produced by the circulating Joscphson currents be- 
tween the grains. As is known p^Jl7j , a constant electric 
field E applied to a single Josephson junction ( JJ) causes 
a time evolution of the phase difference. In this particular 
case Eq.(2) reads 



2e - 
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The resulting AC superconducting current in the junc- 
tion is 



2eEj 



■ sin i 
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If, in addition to the external electric field E, the net- 
work of superconducting grains is under the influence of 
an applied magnetic field H, the frustration parameter 
Aij (t) in Eq. (3) takes the following form 
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Here, Rij = (r*j + r*j)/2, and we have used the conven- 
tional relationship between the vector potential A and a 
constant magnetic field H = rot A (with dH/dt = 0), as 
well as a homogeneous electric field E — —dA/dt (with 
rotE = 0). In the type II HTCS the magnetic self-field 
effects for the array as a whole are expected to be negligi- 
ble . The grains themselves are in fact larger than the 
London penetration depth and we must then have that 
the corresponding Josephson penetration length must be 
much larger than the grain size (since the self-induced 
magnetic fields can in principle be quite pronounced 
for large-size junctions even in zero applied magnetic 
fields |15|| .) Specifically, this is justified for short junc- 
tions with the size d <C Xj, where Xj — \J < f ) o/47r^oj c AL 
is the Josephson penetration length with Xl being the 
grain London penetration depth and j c its Josephson 
critical current density. In particular, since in HTCS 
Xl — 150nm, the above condition will be fulfilled for 
d ~ lfim and j c ~ 10 4 A/m 2 which are the typical pa- 
rameters for HTCS ceramics [^6|. Likewise, to ensure 
the uniformity of the applied electric field, we also as- 
sume that d ^ Xe, where A# is an effective electric field 
penetration depth ||. 

When the AC supercurrent I*j(t) (defined by Eqs.(2), 
(5) and (6)) circulates around a set of grains, that form a 
random area plaquette, it induces a random AC magnetic 
moment jl(t) of the Josephson network jl2| 
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Notice that in the MEE-like effect discussed here for a 
granular superconductor, the electric-field induced mag- 
netic moment in the system is still present in zero ap- 
plied magnetic field due to the phase coherent currents 
between the weakly-coupled superconducting grains. 

To consider the essence of the superconducting elec- 
tric field-induced MEE, we assume that in a zero elec- 
tric field the phase difference between the adjacent grains 
</>ij(0) = which corresponds to a fully coherent state 
of the array. In this particular case, the electric-field in- 
duced averaged magnetization reads 
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M s (E)=p s (t) = - J dt J ,lr,, l dl!,, ! S\r,, r /i 1 , 
o o 
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where r is the electronic relaxation scattering time, and S 
is the joint probability distribution function (see below). 

To obtain an explicit expression for the electric-field 
dependent magnetization, we consider a site positional 
disorder that allows for small random radial displace- 
ments. Namely, the sites in a 3D cubic lattice are as- 
sumed to move from their equilibrium positions accord- 
ing to the normalized (separable) distribution function 



(9) 



It can be shown that the main qualitative results of 
this paper do not depend on the particular choice of 
the probability distribution function. For simplicity here 
we assume an exponential distribution law for the dis- 
tance between grains, P r (r) = P{x\)P{x2)P{xz) with 
P r (xj) — {\/d)e~ Xj / d , and a short range distribution for 
the dependence of the center-of-mass probability Pr(R) 
(around some constant value D). The specific form of 
the latter distribution will not affect the qualitative na- 
ture of the final result. (Notice that in fact the former 
distribution function P r (r) reflects a short-range charac- 
ter of the Josephson coupling in granular superconduc- 
tor. Indeed, according to the conventional picture 
the Josephson coupling Ej(fij) can be assumed to vary 
exponentially with the distance fij between neighboring 
grains, i.e., Ej(fij) — Eje~ K ' rij . For isotropic arrange- 
ment of identical grains, with spacing d between the cen- 
ters of adjacent grains, we have k = (3,3,3) and thus 
d is of the order of an average grain size.) Taking the 

applied electric field along the x-axis, E = (E x ,0, 0), we 
get finally 



M Z (E X ) = ^^--H Z (E X ), 
Mo 



(10) 



for the induced transverse magnetization (along the X3 
z-axis), where 



B Z (E X ) = /iqMq 



E x /E 



1 + (E x /E ) 2 



and 



H Z {E X ) 
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stand for the electric-field induced magnetic induction 
B Z (E X ) and magnetic field H Z {E X ), respectively. The 
induced Josephson current I(E X ) is simply given by 
Ampere's law I(E X ) — H z (E x )d. In these equations, 
Mq = 2neEjNdD /%, with N the total number of grains 
and Eq = H/2der. Eq.(10) is the main result of this 
paper, which we proceed to analyze below. 
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FIG. 1. The induced magnetization m(e) — M z /Mo as a 
function of the normalized applied electric field e = E x /Eo, 
according to Eq.(10). The inset shows a blow up of the para- 
magnetic region of m(e) (exhibiting its linear superconducting 
magnetoelectric effect regime). 

As is seen from Eq.(10), the behavior of the magneti- 
zation in the applied electric field is determined by the 
competition between the two contributions, the magnetic 
induction B Z (E X ) and the current induced magnetic field 
H Z (E X ) (or the corresponding Josephson current I(E X )). 
Namely, below a critical (threshold) field E c m 1.94 Eq 
where B Z (E X ) > iaqH z (E x ), a paramagnetic phase of 
the MEE is due to the modification of the magnetic in- 
duction in the applied electric field. On the other hand 
above the threshold field (when E x becomes larger than 
E c ) the Josephson current I(E X ) induced contribution 
starts to prevail, leading to the appearance of the dia- 
magnetic signal (seen as a small negative part of the in- 
duced magnetization in Fig.l). Such an electric field in- 
duced paramagnetic-to-diamagnetic transition has been 
actually observed [|l]-[| for behavior of the critical cur- 
rent in ceramic HTCS and it was attributed |(| to a 
proximity-mediated enhancement of the superconductiv- 
ity in a granular material in a strong enough electric field. 
Explicitly, the critical current in a YBCO sample was 
found to reach a maximum at E = 4 x l0 7 V/m. To 
relate this experimental value with the model parame- 
ters, first of all, we need to estimate an order of magni- 
tude of the relaxation time t in a zero applied electric 
field. This will provide an upper limit for the relevant 
r-distribution in our system. It is reasonable to connect 
zero-field r = r(0) with the Josephson tunneling time |l(| 
tj = (Ro/Rn)(h/Ej) (where R = h/4e 2 , and R n is the 
normal state resistance between grains). Typically, for 
HTCS ceramics Ej/k B ^ 90K and R n /R ^ 1CT 3 , so 
that tj ~ 10 _10 s. At the same time, at high enough elec- 
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trie fields where the MEE becomes strongly nonlinear, 
we can expect quite a tangible decrease of the relaxation 
time. Indeed, for an average grain size d w lfim, the 
characteristic field Eq — 4 x 10 7 V/m (which corresponds 
to the region where a prominent enhancement of the criti- 
cal current was observed 1|) introduces a substantially 
shorter relaxation time t(Eq) = %/2deEo ~ 10~ 16 s, in 
agreement with observations 

To estimate the relative magnitude of the supercon- 
ducting analog of the MEE predicted here, we can com- 
pare it with the normal (Ohmic) contribution to the mag- 
netization 
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M n (E) = jl n {t) 
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dt / dr i ^dR l ^S{r l j, R hj )p, n , 
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where 



2 « = 7T J2 I Zj 



(13) 
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Here 1^ = V hj /R r , 



is the normal current component due 



to the applied electric field E, with V* -■ = E 



being 



the induced voltage, and R n the normal state resistance 
between grains. As a result, the normal state contribu- 
tion (for E along the rr-axis) reads M n = a n E x , with 
a n = nd 2 DN/R n . Similarly, according to Eq.(lO), the 
low field contribution to the superconducting MEE gives 
M s ~ a s E x with a s = 2ne 2 EjNT{0)d 2 D/h 2 . Thus, at 
low enough applied fields (when E x <£"o) 
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where tj = (R /R n )(h/Ej) with R = h/4e 2 . According 
to our previous discussion on the relevant relaxation-time 
distribution spectrum in our model system, we may con- 
clude that r(0) < tj. So, we arrive at the following ratio 
between the coefficients of the superconducting to normal 
MEEs, namely a s /a n < 1. The above estimate of the 
weak-links induced MEE (along with its rather specific 
field dependence, see Fig.l) suggests quite an optimistic 
possibility to observe the predicted effect experimentally 
in HTCS ceramics or in a specially prepared system of 
arrays of superconducting grains. 

We note that in the present analysis we have not explic- 
itly considered the polarization effects (due to the inter- 
action between the applied electric field and the grain's 
charges) which may become important at high enough 
fields (or for small enough grains), leading to more sub- 
tle phenomena (like Coulomb blockade and reentrant-like 
behavior) that will demand the inclusion of charging en- 
ergy effects in the analysis. 

In summary, we have used a zero-temperature ran- 
dom 3D XY model to predict the appearance of a novel 
electric-field induced magnetization in a granular super- 
conductor (an analog of the magnetoelectric effect). The 



induced magnetization has a very distinctive nonlinear 
and para- to dia-magnetic behavior as the field is in- 
creased. We have also estimated the possible parameter 
range conditions to observe this effect experimentally. 
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